Genetic Diversity Assessment in Plants from Reduced Representation Sequencing Data.

Genetic diversity refers to the variety of genetic traits within a population or a species. It is an essential aspect of both plant ecology and plant breeding because it contributes to the adaptability, survival, and resilience of populations in changing environments. This chapter outlines a pipeline for estimating genetic diversity statistics from reduced representation or whole genome sequencing data. The pipeline involves obtaining DNA sequence reads, mapping the corresponding reads to a reference genome, calling variants from the alignments, and generating an unbiased estimation of nucleotide diversity and divergence between populations. The pipeline is suitable for single-end Illumina reads and can be adjusted for paired-end reads. The resulting pipeline provides a comprehensive approach for aligning and analyzing sequencing data to estimate genetic diversity.

Identification of Novel Induced Mutations in Seed and Vegetatively Propagated Plants from Reduced Representation or Whole Genome Sequencing Data.

Treatment of plants with chemical mutagens results primarily in the production of novel single nucleotide variants. Mutagenesis is a mostly random process and as such plants derived from mutagenesis of different seeds or in vitro material are expected to accumulate different mutations. An important step in the creation of a mutant population for forward or reverse genetics is the choice of treatment conditions (e.g., dosage) such that sufficient mutations accumulate while not adversely affecting propagation of the plant. DNA sequencing provides a quick method to evaluate the effect of different treatment conditions and their effect on the density and spectrum of accumulated mutations. Whole genome sequencing or reduced representation sequencing is carried out followed by mapping to a reference genome and production of a Variant Call Format (VCF) file. We provide here a method for generating a multi-sample VCF from mutagenized plants and describe a new tool to streamline the process of recovering unique induced mutations and determining their possible effect on gene function.

Sample and Library Preparation for PacBio Long-Read Sequencing in Grapevine.

PacBio long-read sequencing is a third-generation technology that generates long reads up to 20 kilobases (kb), unlike short-read sequencing instruments that produce up to 600 bases. Long-read sequencing is particularly advantageous in higher organisms, such as humans and plants, where repetitive regions in the genome are more abundant. The PacBio long-read sequencing uses a single molecule, real-time approach where the SMRT cells contain several zero-mode waveguides (ZMWs). Each ZMW contains a single DNA molecule bound by a DNA polymerase. All ZMWs are flushed with deoxy nucleotides with a fluorophore specific to each nucleotide. As the sequencing proceeds, the detector detects the wavelength of the fluorescence and the nucleotides are read in real-time. This chapter describes the sample and library preparation for PacBio long-read sequencing for grapevine.

Identification of Induced Copy Number Variation from Low Coverage Sequence Data.

Induced mutations have been an important tool for plant breeding and functional genomics for more than 80 years. Novel mutations can be induced by treating seed or other plant cells with chemical mutagens or ionizing radiation. The majority of released mutant crop varieties were developed using ionizing radiation. This has been shown to create a variety of different DNA lesions including large (e.g., >=10,000 bps) copy number variations (CNV). Detection of induced DNA lesions from whole genome sequence data is useful for choosing a mutagen dosage prior to committing resources to develop a large mutant population for forward or reverse-genetic screening. Here I provide a method for detecting large induced CNV from mutant plants that utilizes a new tool to streamline the process of obtaining read coverage directly from BAM files, comparing non-mutagenized controls and mutagenized samples, and plotting the results for visual evaluation. Example data is provided from low coverage sequence data from gamma-irradiated vegetatively propagated triploid banana.

A Comprehensive Protocol for Assembly of Multiple gRNAs into a Direct Vector for Genome Editing in Tomato.

The Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR)-Cas 9 (CRISPR-associated protein 9) is a robust DNA-encoded, RNA-mediated sequence-specific nuclease system widely used for genome editing of various plants. Although there are many reports on the assembly of gRNAs and plant transformation, there is no single resource for the complete gene editing methodology in tomato. This chapter provides a comprehensive protocol for designing gRNAs, their assembly into the vector, plant transformation, and final mutant analysis in tomato.

CRISPR/Cas9-Mediated Genome Editing in Indica Rice (Oryza sativa L. subsp. indica var. CR-5272).

Tissue culture optimization protocols limit indica rice breeding. Such a challenge is vital because emergent techniques still rely on tissue culture methods and could allow the breeding of new varieties with higher production and toleration of adverse environmental effects caused by climate change. Genome editing technology, using CRISPR/Cas9, is a fast and precise method for accelerated plant breeding. It limited its use in indica subspecies because of the recalcitrant response to in vitro culture methods. This chapter describes a protocol for CRISPR/Cas9 editing in indica subspecies, specifically in the CR-5272 variety derived from parental lines IR-822, using Agrobacterium tumefaciens and biolistic transformation.

Genome Editing in Brassica juncea Using CRISPR/Cas9 Technology.

Modern genome editing tools particularly CRISPR/Cas9 have revolutionized plant genome manipulation for engineering resilience against changing climatic conditions, disease infestation, as well as functional genomic studies. CRISPR-mediated genome editing allows for editing at a single as well as multiple locations in the genome simultaneously, making it an effective tool for polyploid species too. However, still, its applications are limited to the model crops only. Extending it to crop plants will help improve field crops against the changing climates more rapidly and precisely. Here we describe the protocol for editing the genome of a field crop Brassica juncea (mustard), an allotetraploid and important oilseed crop of the Indo-Pak Subcontinent region. This protocol is based on the Agrobacterium-mediated transformation for the delivery of CRISPR components into the plant genome using cotyledon as explants. We elaborate on steps for recovering genome-edited knockouts, for validation of the edits, as well as recovering the transgene-free edited plants through a commonly used segregating approach.

Subgenome-aware analyses reveal the genomic consequences of ancient allopolyploid hybridizations throughout the cotton family.

Malvaceae comprise some 4,225 species in 243 genera and nine subfamilies and include economically important species, such as cacao, cotton, durian, and jute, with cotton an important model system for studying the domestication of polyploids. Here, we use chromosome-level genome assemblies from representatives of five or six subfamilies (depending on the placement of Ochroma) to differentiate coexisting subgenomes and their evolution during the family's deep history. The results reveal that the allohexaploid Helicteroideae partially derive from an allotetraploid Sterculioideae and also form a component of the allodecaploid Bombacoideae and Malvoideae. The ancestral Malvaceae karyotype consists of 11 protochromosomes. Four subfamilies share a unique reciprocal chromosome translocation, and two other subfamilies share a chromosome fusion. DNA alignments of single-copy nuclear genes do not yield the same relationships as inferred from chromosome structural traits, probably because of genes originating from different ancestral subgenomes. These results illustrate how chromosome-structural data can unravel the evolutionary history of groups with ancient hybrid genomes.

Deep learning the cis-regulatory code for gene expression in selected model plants.

Elucidating the relationship between non-coding regulatory element sequences and gene expression is crucial for understanding gene regulation and genetic variation. We explored this link with the training of interpretable deep learning models predicting gene expression profiles from gene flanking regions of the plant species Arabidopsis thaliana, Solanum lycopersicum, Sorghum bicolor, and Zea mays. With over 80% accuracy, our models enabled predictive feature selection, highlighting e.g. the significant role of UTR regions in determining gene expression levels. The models demonstrated remarkable cross-species performance, effectively identifying both conserved and species-specific regulatory sequence features and their predictive power for gene expression. We illustrated the application of our approach by revealing causal links between genetic variation and gene expression changes across fourteen tomato genomes. Lastly, our models efficiently predicted genotype-specific expression of key functional gene groups, exemplified by underscoring known phenotypic and metabolic differences between Solanum lycopersicum and its wild, drought-resistant relative, Solanum pennellii.

Genome-wide identification, characterization and expression of C2H2 zinc finger gene family in Opisthopappus species under salt stress.

BACKGROUND: The C2H2 zinc finger protein family plays important roles in plants. However, precisely how C2H2s function in Opisthopappus (Opisthopappus taihangensis and Opisthopappus longilobus) remains unclear. RESULTS: In this study, a total of 69 OpC2H2 zinc finger protein genes were identified and clustered into five Groups. Seven tandem and ten fragment repeats were found in OpC2H2s, which underwent robust purifying selection. Of the identified motifs, motif 1 was present in all OpC2H2s and conserved at important binding sites. Most OpC2H2s possessed few introns and exons that could rapidly activate and react when faced with stress. The OpC2H2 promoter sequences mainly contained diverse regulatory elements, such as ARE, ABRE, and LTR. Under salt stress, two up-regulated OpC2H2s (OpC2H2-1 and OpC2H2-14) genes and one down-regulated OpC2H2 gene (OpC2H2-7) might serve as key transcription factors through the ABA and JA signaling pathways to regulate the growth and development of Opisthopappus species. CONCLUSION: The above results not only help to understand the function of C2H2 gene family but also drive progress in genetic improvement for the salt tolerance of Opisthopappus species.

Telomere-to-telomere genome assembly of Oldenlandia diffusa.

We report the complete telomere-to-telomere genome assembly of Oldenlandia diffusa which renowned in traditional Chinese medicine, comprising 16 chromosomes and spanning 499.7 Mb. The assembly showcases 28 telomeres and minimal gaps, with a total of only five. Repeat sequences constitute 46.41% of the genome, and 49,701 potential protein-coding genes have been predicted. Compared with O. corymbosa, O. diffusa exhibits chromosome duplication and fusion events, diverging 20.34 million years ago. Additionally, a total of 11 clusters of terpene synthase have been identified. The comprehensive genome sequence, gene catalog, and terpene synthase clusters of O. diffusa detailed in this study will significantly contribute to advancing research in this species' genetic, genomic, and pharmacological aspects.

The genome and population genomics of allopolyploid Coffea arabica reveal the diversification history of modern coffee cultivars.

Coffea arabica, an allotetraploid hybrid of Coffea eugenioides and Coffea canephora, is the source of approximately 60% of coffee products worldwide, and its cultivated accessions have undergone several population bottlenecks. We present chromosome-level assemblies of a di-haploid C. arabica accession and modern representatives of its diploid progenitors, C. eugenioides and C. canephora. The three species exhibit largely conserved genome structures between diploid parents and descendant subgenomes, with no obvious global subgenome dominance. We find evidence for a founding polyploidy event 350,000-610,000 years ago, followed by several pre-domestication bottlenecks, resulting in narrow genetic variation. A split between wild accessions and cultivar progenitors occurred ~30.5 thousand years ago, followed by a period of migration between the two populations. Analysis of modern varieties, including lines historically introgressed with C. canephora, highlights their breeding histories and loci that may contribute to pathogen resistance, laying the groundwork for future genomics-based breeding of C. arabica.

Pangenome analysis reveals transposon-driven genome evolution in cotton.

BACKGROUND: Transposable elements (TEs) have a profound influence on the trajectory of plant evolution, driving genome expansion and catalyzing phenotypic diversification. The pangenome, a comprehensive genetic pool encompassing all variations within a species, serves as an invaluable tool, unaffected by the confounding factors of intraspecific diversity. This allows for a more nuanced exploration of plant TE evolution. RESULTS: Here, we constructed a pangenome for diploid A-genome cotton using 344 accessions from representative geographical regions, including 223 from China as the main component. We found 511 Mb of non-reference sequences (NRSs) and revealed the presence of 5479 previously undiscovered protein-coding genes. Our comprehensive approach enabled us to decipher the genetic underpinnings of the distinct geographic distributions of cotton. Notably, we identified 3301 presence-absence variations (PAVs) that are closely tied to gene expression patterns within the pangenome, among which 2342 novel expression quantitative trait loci (eQTLs) were found residing in NRSs. Our investigation also unveiled contrasting patterns of transposon proliferation between diploid and tetraploid cotton, with long terminal repeat (LTR) retrotransposons exhibiting a synchronized surge in polyploids. Furthermore, the invasion of LTR retrotransposons from the A subgenome to the D subgenome triggered a substantial expansion of the latter following polyploidization. In addition, we found that TE insertions were responsible for the loss of 36.2% of species-specific genes, as well as the generation of entirely new species-specific genes. CONCLUSIONS: Our pangenome analyses provide new insights into cotton genomics and subgenome dynamics after polyploidization and demonstrate the power of pangenome approaches for elucidating transposon impacts and genome evolution.

PPanG: a precision pangenome browser enabling nucleotide-level analysis of genomic variations in individual genomes and their graph-based pangenome.

Graph-based pangenome is gaining more popularity than linear pangenome because it stores more comprehensive information of variations. However, traditional linear genome browser has its own advantages, especially the tremendous resources accumulated historically. With the fast-growing number of individual genomes and their annotations available, the demand for a genome browser to visualize genome annotation for many individuals together with a graph-based pangenome is getting higher and higher. Here we report a new pangenome browser PPanG, a precise pangenome browser enabling nucleotide-level comparison of individual genome annotations together with a graph-based pangenome. Nine rice genomes with annotations were provided by default as potential references, and any individual genome can be selected as the reference. Our pangenome browser provides unprecedented insights on genome variations at different levels from base to gene, and reveals how the structures of a gene could differ for individuals. PPanG can be applied to any species with multiple individual genomes available and it is available at https://cgm.sjtu.edu.cn/PPanG .

Long-read sequencing of extrachromosomal circular DNA and genome assembly of a Solanum lycopersicum breeding line revealed active LTR retrotransposons originating from S. Peruvianum L. introgressions.

Transposable elements (TEs) are a major force in the evolution of plant genomes. Differences in the transposition activities and landscapes of TEs can vary substantially, even in closely related species. Interspecific hybridization, a widely employed technique in tomato breeding, results in the creation of novel combinations of TEs from distinct species. The implications of this process for TE transposition activity have not been studied in modern cultivars. In this study, we used nanopore sequencing of extrachromosomal circular DNA (eccDNA) and identified two highly active Ty1/Copia LTR retrotransposon families of tomato (Solanum lycopersicum), called Salsa and Ketchup. Elements of these families produce thousands of eccDNAs under controlled conditions and epigenetic stress. EccDNA sequence analysis revealed that the major parts of eccDNA produced by Ketchup and Salsa exhibited low similarity to the S. lycopersicum genomic sequence. To trace the origin of these TEs, whole-genome nanopore sequencing and de novo genome assembly were performed. We found that these TEs occurred in a tomato breeding line via interspecific introgression from S. peruvianum. Our findings collectively show that interspecific introgressions can contribute to both genetic and phenotypic diversity not only by introducing novel genetic variants, but also by importing active transposable elements from other species.

Genetic insights into agronomic and morphological traits of drug-type cannabis revealed by genome-wide association studies.

Cannabis sativa L., previously concealed by prohibition, is now a versatile and promising plant, thanks to recent legalization, opening doors for medical research and industry growth. However, years of prohibition have left the Cannabis research community lagging behind in understanding Cannabis genetics and trait inheritance compared to other major crops. To address this gap, we conducted a comprehensive genome-wide association study (GWAS) of nine key agronomic and morphological traits, using a panel of 176 drug-type Cannabis accessions from the Canadian legal market. Utilizing high-density genotyping-by-sequencing (HD-GBS), we successfully generated dense genotyping data in Cannabis, resulting in a catalog of 800 K genetic variants, of which 282 K common variants were retained for GWAS analysis. Through GWAS analysis, we identified 18 markers significantly associated with agronomic and morphological traits. Several identified markers exert a substantial phenotypic impact, guided us to putative candidate genes that reside in high linkage-disequilibrium (LD) with the markers. These findings lay a solid foundation for an innovative cannabis research, leveraging genetic markers to inform breeding programs aimed at meeting diverse needs in the industry.

Analyzing genetic diversity in luffa and developing a Fusarium wilt-susceptible linked SNP marker through a single plant genome-wide association (sp-GWAS) study.

BACKGROUND: Luffa (Luffa spp.) is an economically important crop of the Cucurbitaceae family, commonly known as sponge gourd or vegetable gourd. It is an annual cross-pollinated crop primarily found in the subtropical and tropical regions of Asia, Australia, Africa, and the Americas. Luffa serves not only as a vegetable but also exhibits medicinal properties, including anti-inflammatory, antidiabetic, and anticancer effects. Moreover, the fiber derived from luffa finds extensive applications in various fields such as biotechnology and construction. However, luffa Fusarium wilt poses a severe threat to its production, and existing control methods have proven ineffective in terms of cost-effectiveness and environmental considerations. Therefore, there is an urgent need to develop luffa varieties resistant to Fusarium wilt. Single-plant GWAS (sp-GWAS) has been demonstrated as a promising tool for the rapid and efficient identification of quantitative trait loci (QTLs) associated with target traits, as well as closely linked molecular markers. RESULTS: In this study, a collection of 97 individuals from 73 luffa accessions including two major luffa species underwent single-plant GWAS to investigate luffa Fusarium wilt resistance. Utilizing the double digest restriction site associated DNA (ddRAD) method, a total of 8,919 high-quality single nucleotide polymorphisms (SNPs) were identified. The analysis revealed the potential for Fusarium wilt resistance in accessions from both luffa species. There are 6 QTLs identified from 3 traits, including the area under the disease progress curve (AUDPC), a putative disease-resistant QTL, was identified on the second chromosome of luffa. Within the region of linkage disequilibrium, a candidate gene homologous to LOC111009722, which encodes peroxidase 40 and is associated with disease resistance in Cucumis melo, was identified. Furthermore, to validate the applicability of the marker associated with resistance from sp-GWAS, an additional set of 21 individual luffa plants were tested, exhibiting 93.75% accuracy in detecting susceptible of luffa species L. aegyptiaca Mill. CONCLUSION: In summary, these findings give a hint of genome position that may contribute to luffa wild resistance to Fusarium and can be utilized in the future luffa wilt resistant breeding programs aimed at developing wilt-resistant varieties by using the susceptible-linked SNP marker.

Dynamic changes in the plastid and mitochondrial genomes of the angiosperm Corydalis pauciovulata (Papaveraceae).

BACKGROUND: Corydalis DC., the largest genus in the family Papaveraceae, comprises > 465 species. Complete plastid genomes (plastomes) of Corydalis show evolutionary changes, including syntenic arrangements, gene losses and duplications, and IR boundary shifts. However, little is known about the evolution of the mitochondrial genome (mitogenome) in Corydalis. Both the organelle genomes and transcriptomes are needed to better understand the relationships between the patterns of evolution in mitochondrial and plastid genomes. RESULTS: We obtained complete plastid and mitochondrial genomes from Corydalis pauciovulata using a hybrid assembly of Illumina and Oxford Nanopore Technologies reads to assess the evolutionary parallels between the organelle genomes. The mitogenome and plastome of C. pauciovulata had sizes of 675,483 bp and 185,814 bp, respectively. Three ancestral gene clusters were missing from the mitogenome, and expanded IR (46,060 bp) and miniaturized SSC (202 bp) regions were identified in the plastome. The mitogenome and plastome of C. pauciovulata contained 41 and 67 protein-coding genes, respectively; the loss of genes was a plastid-specific event. We also generated a draft genome and transcriptome for C. pauciovulata. A combination of genomic and transcriptomic data supported the functional replacement of acetyl-CoA carboxylase subunit ß (accD) by intracellular transfer to the nucleus in C. pauciovulata. In contrast, our analyses suggested a concurrent loss of the NADH-plastoquinone oxidoreductase (ndh) complex in both the nuclear and plastid genomes. Finally, we performed genomic and transcriptomic analyses to characterize DNA replication, recombination, and repair (DNA-RRR) genes in C. pauciovulata as well as the transcriptomes of Liriodendron tulipifera and Nelumbo nuicifera. We obtained 25 DNA-RRR genes and identified their structure in C. pauciovulata. Pairwise comparisons of nonsynonymous (dN) and synonymous (dS) substitution rates revealed that several DNA-RRR genes in C. pauciovulata have higher dN and dS values than those in N. nuicifera. CONCLUSIONS: The C. pauciovulata genomic data generated here provide a valuable resource for understanding the evolution of Corydalis organelle genomes. The first mitogenome of Papaveraceae provides an example that can be explored by other researchers sequencing the mitogenomes of related plants. Our results also provide fundamental information about DNA-RRR genes in Corydalis and their related rate variation, which elucidates the relationships between DNA-RRR genes and organelle genome stability.

Mining genomic regions associated with agronomic and biochemical traits in quinoa through GWAS.

Quinoa (Chenopodium quinoa Willd.), an Andean crop, is a facultative halophyte food crop recognized globally for its high nutritional value and plasticity to adapt to harsh conditions. We conducted a genome-wide association study on a diverse set of quinoa germplasm accessions. These accessions were evaluated for the following agronomic and biochemical traits: days to 50% flowering (DTF), plant height (PH), panicle length (PL), stem diameter (SD), seed yield (SY), grain diameter (GD), and thousand-grain weight (TGW). These accessions underwent genotyping-by-sequencing using the DNBSeq-G400R platform. Among all evaluated traits, TGW represented maximum broad-sense heritability. Our study revealed average SNP density of ≈ 3.11 SNPs/10 kb for the whole genome, with the lowest and highest on chromosomes Cq1B and Cq9A, respectively. Principal component analysis clustered the quinoa population in three main clusters, one clearly representing lowland Chilean accessions, whereas the other two groups corresponded to germplasm from the highlands of Peru and Bolivia. In our germplasm set, we estimated linkage disequilibrium decay to be ≈ 118.5 kb. Marker-trait analyses revealed major and consistent effect associations for DTF on chromosomes 3A, 4B, 5B, 6A, 7A, 7B and 8B, with phenotypic variance explained (PVE) as high as 19.15%. Nine associations across eight chromosomes were also found for saponin content with 20% PVE by qSPN5A.1. More QTLs were identified for PL and TGW on multiple chromosomal locations. We identified putative candidate genes in the genomic regions associated with DTF and saponin content. The consistent and major-effect genomic associations can be used in fast-tracking quinoa breeding for wider adaptation across marginal environments.

Genome-wide characterization of DNA methyltransferase family genes implies GhDMT6 improving tolerance of salt and drought on cotton.

BACKGROUND: DNA methylation is an important epigenetic mode of genomic DNA modification and plays a vital role in maintaining epigenetic content and regulating gene expression. Cytosine-5 DNA methyltransferase (C5-MTase) are the key enzymes in the process of DNA methylation. However, there is no systematic analysis of the C5-MTase in cotton so far, and the function of DNMT2 genes has not been studied. METHODS: In this study, the whole genome of cotton C5-MTase coding genes was identified and analyzed using a bioinformatics method based on information from the cotton genome, and the function of GhDMT6 was further validated by VIGS experiments and subcellular localization analysis. RESULTS: 33 C5-MTases were identified from three cotton genomes, and were divided into four subfamilies by systematic evolutionary analysis. After the protein domain alignment of C5-MTases in cotton, 6 highly conserved motifs were found in the C-terminus of 33 proteins involved in methylation modification, which indicated that C5-MTases had a basic catalytic methylation function. These proteins were divided into four classes based on the N-terminal difference, of which DNMT2 lacks the N-terminal regulatory domain. The expression of C5-MTases in different parts of cotton was different under different stress treatments, which indicated the functional diversity of cotton C5-MTase gene family. Among the C5-MTases, the GhDMT6 had a obvious up-regulated expression. After silencing GhDMT6 with VIGS, the phenotype of cotton seedlings under different stress treatments showed a significant difference. Compared with cotton seedlings that did not silence GhDMT6, cotton seedlings silencing GhDMT6 showed significant stress resistance. CONCLUSION: The results show that C5-MTases plays an important role in cotton stress response, which is beneficial to further explore the function of DNMT2 subfamily genes.